First principle computation of stripes in cuprates 
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We present a first principle computation of vertical stripes in Lai5/gSri/gCu04 witliin tlie LDA+J7 
metfiod. We find tliat Cu centered stripes are unstable toward O centered stripes. The metallic core 
of the stripe is quite wide and shows reduced magnetic moments with suppressed antiferromagnetic 
(AF) interactions. The system can be pictured as alternating metallic and AF two-leg ladders the 
latter with strong AF interaction and a large spin gap. The Fermi surface shows warping due to 
interstripe hybridization. The periodicity and amplitude of the warping is in good agreement with 
angle resolved photoemission experiment. We discuss the connection with low-energy theories of 
the cuprates. 
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Recent theories of the high temperature superconduct- 
ing cuprates are based on the idea that hole-rich quasi 
one-dimensional textures (stripes) are the basic buildin 
blocs to understand the low-energy physics 
ISIllilSIISIIlIIlIIi- These model computations ei- 
ther assume stripes to obtain a low-energy quantum field 
theory (QFT)|1 H H H S [H JTior obtain stripes 



in simplified model Hamiltoniansj^TS 1^ Qi 

For the QFT's detailed information about the symme- 
try and the extent of real stripes is lacking so much 
of the modeling is based on guesses and the relevant 
regimes are identified a posteriori, by confrontation with 
experiment 0, 0, 0, 0, leading to a plethora of 
possible scenarios. 

In order to constrain low-energy models some progress 
can be made by solving more or less realistic Hamiltoni- 
ans in various approximations |Q 0, 0, 0, ll^, 0, 113, IU 
however these are tight to strong assumptions about 
the electronic structure since are restricted to a small 
number of orbitals and/or short range interactions and 
a single Cu-0 plane. Here we present a first principle 
computation of stripes in cuprates based on the LDA-I-J7 
method for doping x = 1/8 and the experimentally mea- 
sured magnetic incomensurability e = 1/8. Our compu- 
tation are based on the full orbital variational space and 
the three-dimensional Coulomb interaction and so go be- 
yond the previous works. LDA based methods are widely 
believed to give reliable charge distributions and there- 
fore can provide a realistic picture of stripe textures. We 
find that a Cu centered solution is not stable. Instead the 
stripe is O centered and consist of two columns of "metal- 
lic" Cu (the core) and two columns of antiferromagnetic 
(AF) Cu forming an alternating system of metallic and 
AF two-leg ladders. 

Added holes are spread over the "metallic" Cu rows 
and the surrounding O leading to a stripe core much 
wider than the traditional stripe picture [l^. The 



electronic structure consists of a half-filed quasi one- 
dimensional band for each stripe core corresponding to 
motion along the stripe. The system can be schematized 
as an array of one-dimensional electron gases (IDEG) 
separated by spin-ladders ( "the environment" in the lan- 
guage of Ref. 0, 0, S 0) with a large spin-gap. Inter- 
stripe hybridization is not negligible but leads to a notice- 
able warping of the Fermi surface which is in good agree- 
ment with angle resolved photoemission experiments[25|. 
Additionally there is approximate particle-hole symme- 
try around the Fermi level in good agreement with trans- 
port experiments H^l- We provide estimates of the 
relevant energy scales for "striped" low-energy theories 
of cuprates. 
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FIG. 1: The hole-charge and spin order within CUO2 plane. 
Sites with (without) arrows (indicating spin directions) are Cu 
(O). The radii of circles are proportional to the hole content 
inside atomic spheres. The grey regions are metallic core with 
Cu of type Cul and the white regions are the AF ladders with 
Cu of type Cu2. Numbers indicate the distribution of the 
doped hole. Not shown are 1,2 % (0,6 %) of hole charge in 
the apical O's belonging to Cul (Cu2) and 0,3 % of charge on 
La. The unit cell in the CUO2 plane is shown by the dashed 
line. 

In the present work we used a LDA-t-fJ approach |25 
Egf realized in frame of scalar relativistic LMTO 
method within the Atomic Sphere Approximation 
(ASA) [13 with the Coulomb term of LDA-hC/ func- 
tional proposed in Ref. |2^. The radii of atomic 
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FIG. 2: Wannier function of the active band. 



spheres were r(Cu)=2.2 a.u., r(0-plain)=1.94 a.u., r(0- 
apex)=2.00 a.u., r(La)=3.2 a.u; three types of empty 
spheres were added. The values of on-cite Coulomb U 
and exchange J parameters were chosen 8 and 0.88 eV, 
correspondingly. 

We studied a supercell consisted of eight formula units 
with an additional hole in the absence of the Sr impu- 
rity, compensating the charge by adding suitable nega- 
tive charges to all atoms in the supercell. A Cu centered 
solution could not be stabilized, not even as a metastable 
state. In Fig. ^ we show the charge and spin distribution 
of the stable oxygen centered solution. The hole charge 
is concentrated on the two core Cu's forming the AF 
domain wall and the surrounding oxygens (grey region) 
and leak significantly over the AF regions indicating that 
interstripe overlap is not negligible as discussed below. 

The stripe core consist of rungs of ferromagnetic Cu 
sites. This locally enhances the effective Cu-Cu hy- 
bridization and splits states into the charge transfer band 
which disperse strongly in the direction of the stripe 
(Fig. I2J). The band crossing the Fermi level (here- 
after "the active band"), is half- filled, has quasi one- 
dimensional character and is similar to the one found in 
the three-band Hubbard model|0. The electronic struc- 
ture is approximately particle-hole symmetric around the 
chemical potential. This is in agreement with the pro- 
posal that approximate particle-hole symmetry in the 
stripe state is the reason why Hall effect [iil and the 
thermoDOw erl2fj t end to be suppressed in the presence 

of stripes[i[3l2lll3- 

A popular model to study the stability and sym- 
metry of stripes in cuprates is the single-band Hub- 
bard model and its strong coupling version, the t — J 
modelUllillii. These models are usually derived by 
considering the low-energy many-body states of a siiigle 
hole added to the insulator. Zhang and Rice (ZR)|35j 
pointed out that the hole will occupy a particular lin- 
ear superposition of the four p oxygen orbitals pointing 
toward a Cu that has the same symmetry as the Cu 
dx2_y2 orbital. The hole in this "ZR-0 orbital" and the 
hole on the Cu form a ZR singlet. The other (orthogo- 



nal) combinations of p oxygen orbitals are projected out. 
By doing so the ZR-0 orbital is not allowed to polar- 
ize. Such approximation is justified in the dilute limit 
however when one considers the stability of dense phases 
special care must be taken on defining low-energy mod- 
els. In particular for CDW phases, polarization involving 
high energy states can significantly affect the stability 
and symmetry Indeed this problem involves the 
short range (high energy) physics of the system and the 
above procedure may be inadequate. 

To illustrate the problem we have computed the Wan- 
nier function (WF) of the active band (Fig.EJ. The WF 
is centered on the O forming the stripe core but extends 
considerably to the surrounding oxygens. It is clear that 
the Cu-centered ZR-0 orbitals are strongly polarized and 
so the usual Hubbard or t — J mapping becomes unre- 
liable. This is in contrast with what was found for one 
hole added to the AF in LDA+C/ or inhomogeneous HF 
where the solution did involve a ZR-0 orbital Is^ . 
We conclude that the stability of stripes can not be re- 
liably studied within the t — J or single-band Hubbard 
model. (Of course the possibility of stripes in those mod- 
els is a problem of interests in itself but does not have 
direct experimental relevance). In particular configura- 
tions centered in the O will be unfavored since involve a 
strong O polarization and single band models neglect the 
relaxation energy involved. If one refrains from studding 
the stability of stripes one can still assume stripes and 
use an effective model to describe the low-energy physics. 

Coulomb effects (including the polarization effects 
mention above) produce a splitting of the Madelung po- 
tential for a hole in a Cu in the metallic region and a Cu in 
the AF region of 0.2eV favoring the charge segregation. 
In LDA this splitting is static but a truly many-body 
treatment should take into account the dynamical effects 
of the polarization. 

Other parameters are similar to the ones in the dilute 
limit |40l |41|. i.e. The metallic regions can be described 
by a two-leg Hubbard model with a hopping along the 
legs and along the rungs t\\ = t± = 0.52 eV and a hopping 
along the diagonals t' = —0.08 eV. By treating this model 
in mean-field with a filling of 3/4 electrons and fitting 
the gap to the upper bands we obtain U = 4.5 eV for the 
Hubbard on-site interaction, similar to the value in the 
dilute case (~ 5 eV) 40J. 

The active band is formed by the odd combination of 
the leg orbitals. The associated Wannier functions are 
quite different for odd and even combinations, in partic- 
ular only the odd combination can have admixture with 
the central px orbital as shown in Fig. |21 which shows 
again that in terms of ZR-0 orbitals large relaxations 
are involved not taken into account in the standard Hub- 
bard model which weights all O's around a Cu on an 
equal manner. 

The two bands immediately below the active one are 
formed by the orbitals of the AF ladders and are com- 



3 




FIG. 3: Calculated band structure of Lai5/gSri/8Cu04 
around the Fermi level {Ep=0). Left (right) panel corre- 
sponds to the direction parallel (perpendicular) to the stripe. 
The point B corresponds to momentum [0, n/{2a)] in the Cu- 
O plane. 



pletely filled in agreement with the assumption that the 
environment is an insulator 0, 0, 0, 0- Therefore the 
AF regions can be viewed as AF ladders. We computed 
the magnetic interaction parameters according to pre- 
scription of Ref. ji^. The interaction along the rungs, 
Jj_ = 146 meV, results to be enhanced respect to the 
value obtained by doing the same computation in the 
stoichiometric system J = 109 mevClj whereas the in- 
teraction along the legs, J|| = 78 meV, is slightly reduced. 
If one consider the leg isolated from the IDEG this would 
lead to a large spin gap of order ~ 95 meV[43i|. 

Coming back to the metallic regions, it is believed that 
transverse hopping will drive an array of IDEG to a two 
dimensional Fermi liquid or an ordered state 45j . How- 
ever according to Ref. H, H, Hi when one couples the 
magnetic system and the IDEG a spin gap will be in- 
duced in the IDEG. This effect will block the effective 
transverse hopping but not the pair hopping0,0,0|. The 
system can become a sliding Luttinger phase and even- 
tually a superconductor. Other related scenarios include 
a boson-fermion model where the bosons are composed 
fermions coupled antiferromagnetically in the insulating 
regions'l4] and a phase with bond order but with two- 
dimensional fermions ,9,„ .lOj. The large values of the AF 
interaction in the insulating regions favor all these sce- 
narios. However to make more progress an estimate of 
the hybridization perpendicular to the stripes is essential. 

The crossing of the active band in the F-B direction 
in Fig. 13 does not occur exactly at the middle point 
(as would be the case for a perfectly half-filled one- 
dimensional band) but is shifted to the B point. This 
shift depends on the momentum perpendicular to the 
stripe and is precisely due to interstripe hopping, i.e 
holes from the core of one stripe can tunnel to the next 
stripe through the AF insulating like region. As a re- 
sult the Fermi surface is not flat but it warps as shown 
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FIG. 4: Angle resolved photoemission spectral weight inte- 
grated within 500 meV of the Fermi level, as a function of 
kx and fc„(25>] together with calculated Fermi surface for 2D 
Cu02 plane with stripes oriented along y (dashed line) . 
Notice that the experimental Fermi surface has contributions 
due to stripes oriented along y and along x the latter produc- 
ing the vertical structures at kx ~ 0.27r/a. 



in Fig. ^ Interestingly this explains oscillations exhibit 
on the experimental Fermi surface in the presence of 
stripeS|25i|. In Fig. 0] we show the experimental and the 
theoretical data. The amplitude and the periodicity of 
the wrapping is in good agreement with experiment. The 
large Fermi surface encounter at larger doping has been 
shown to be in agreement with LDAji^l- The present 
result shows that in underdoped materials LDA based 
methods can also predict the Fermi surface including fine 
details. 

The Fermi surface wrapping can be modeled by adding 
an effective hybridization t_L,cff from core to core to the ef- 
fective two-leg Hubbard model mention above. This gives 
a contribution to the single electron dispersion of the 
form — 2i^.off cos(4afca;) with t±,cS = 0.015 eV. Notice 
that this effect does not break the approximate particle- 
hole symmetry of the system relevant for transport ex- 
periments. The effective hybridization at the Fermi sur- 
face will be of course renormalized by many-body effects 
of the kind considered in Ref. H, |1{ which will tend 
to suppress it or hopping with spin- wave fluctuations |46| 
which will tend to enhance it. The fact that the wrap- 
ping observed is close to the experimental one shows that 
our estimate is not far from the final value although one 
should take into account that the experiment is not done 
strictly at zero energy but integrating in a finite window. 
More theoretical work is needed to decide the nature of 
the ground state with this value of the transverse hop- 
ping. 

To conclude we have presented an LDA+U study of 
vertical stripe in cuprates. Stripes result to be bond- 
centered and quite wide in contrast with the traditional 
stripe picture [ig. The electronic structure shows ap- 
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proximate particle-hole symmetry around the Fermi level 
in good agreement with the picture deduced from trans- 
port experiments [2^ E^ . The system can be pictured 
as alternating metallic Hubbard and insulating Heisen- 
berg ladders. O polarization effects are important and, 
at the static level, they contribute to a large Madelung 
potential splitting between insulating and metallic re- 
gions. The dynamical consequences should be further 
explored. A large exchange interaction in the insulating 
regions favor mechanisms for superconductivity based on 
preformed AF singlets which can induce superconductiv- 
ity on the metallic regions. A non-negligibly hopping 
transverse from stripe to stripe exist which produces ob- 
servable effects and is at the core of the low-energy QFT 
theories. This produces wrapping of the Fermi surface 
which is in agreement with experimental observation. We 
believe numerical studies of simplified models taking into 
account these ingredients should shed light on the prob- 
lem of superconductivity. 
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neighboring Cu02 layers whereas experimentally they was niinirnized by taking a cut of the three-dimensional 

are staked perpendicularly. This produced an unreliable Fermi-surface in a plane where the c-hybridization van- 

c axis hibrydization. The effect on the Fermi surface ish. 



